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Abstract 

We obtain an asymptotic evaluation of the integral 

pCO 

/ e~^ dx 

J -^ 271 - 1-1 

for n —> 00 , where Hn{x) is the Hermite polynomial. This integral is used to deter¬ 
mine the probability for the quantum harmonic oscillator in the nth energy eigenstate 
to tunnel into the classically forbidden region. Numerical results are given to illus¬ 
trate the accuracy of the expansion. 
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1. Introduction 

The quantum harmonic oscillator in dimensionless variables has the Hamiltonian given 
hy H = + x"^), where x,p are the position and momentum operators satisfying 

[x,p\ = i. The normalised eigenstates are 

where Hn{x) is the nth Hermite polynomial. The corresponding probability densities 
Pn{x) are then given by |V’n(x)p, that is 

Pn{x) = ane-^^Hl{x), an = 

The classical turning points are situated at x = ±\/2n -|- 1, so that the probability of 
tunnelling into the classically forbidden region is expressed as 

j-OO 2 

Pn,tun — 2anQni Qn — / _6 H^{x)dx. (f-^) 

Jy/2n+l 
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It has been pointed out recently by Jadczyk [1] that the tunnelling probabilities Pn,tun 
are rarely discussed in the literature on quantum mechanics. In [1], he has shown that 
Pn,tun has the approximate representation for large values of n given by 

Pn,tun = ^ |o.l33975 - + 0(n-^)| (n ^ oo). (1.3) 

Our aim in this note is to obtain more precise information on the large-n behaviour of 
Pn,tun- We carry this out by exploiting the uniform asymptotic expansion of the parabolic 
cylinder function U{—n — 2^x) given in [4], valid for n —)■ oo and x > \/2n + 1 . 


2. A representation of Qn 

We set V := ^/2n + 1 and rescale the variable x to write the integral Qn in (1.2) as 

Qn = ^ J H^{ux) dx. 

The integrand can be expressed in terms of the parabolic cylinder function U{a, z) by [3, 
Eq. (18.15.28)] 

e-(''^)V^Hn{px) = 2^’'"-^U{-^u‘^,2hx). 

The asymptotic expansion of this last function in terms of the Airy function Ai( 2 ;) 
and its derivative is given by [4, Eq. (3.23)] 




U{-hu ,22 ux) = 


il22 + i 
1/2 ^6 


x^ — 1 


T.(C), 


T,(C) = Ai F,{C) + G,(C) 

as z/ —7- +00 uniformly valid for x > 1, where 


^3/2 _ — 1 — I arccoshx (x > 1). 

The functions Fu{C) and Gy{Q) have the asymptotic expansions 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


FuiO ~ E G.iO ~ E 9k{C> 

k=0 k=0 


-2k 


where, from [4, §3.2], we have 


/o(C) — 1) /i(C) — / 2 (C) — ai(C) + 


go{0 = boiC), 9i{0 = 


boiC) 

24 ^ 


with 


1 f 145 + 249x2 - 9x^ 7x(x2 - 6) 455 ) 

“ IIM I (x2 - 1)3 (x2 _ l)3/2^3/2 " 4 ^/’ 


(2.4) 
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2(5 1 12(^2 - 1)3/2 ' 24(3/2 J • 


(2.5) 


Then, from (2.1) and (2.2) it follows that 

22 n+l(^!) 2 gn+i .00 / ^ 


Qn — 


r'2"'+3 


x2 — 1 


1/2 


TliOdx. 


The quantity ( is a monotonically increasing function for x > 1 and ( = 0 when x = 1. 
Furthermore 

dC_( ( 

dx \x2 — 1 

so that we can introduce the new integration variable ( to find 


^ 22-+i(n!)2e"+i 

Qn — / 4^\C)^u{C) dCj 


m ■■= 


( 


X2 — 1 


( 2 . 6 ) 


From (2.2), we have 

T2(() = Ai2(z.4/3^) ir2(^) + 2u-8/3Ai (u^/3()Ai'(u"/3() F,(()G,(() 

1 + 1152/2(0 


= AiV+)(l + ^ + 


5761/4 


+ o(u-o) 


+u-8/36o(()[Ai2 (0/3()]'(^l + + O(u-O) + i/-''*/36g(()Ai'2(0/3()(l + 0(z/-2)). 

Retaining terms up to and including z/“44/3^ therefore find 

J^<PiC)'Tl{C)dC = |^'"(.((){AiOu"/\)(l + Y^+ 

(2.7) 


576^4 J 
+iy-^/\{C) [Ai2 (0/3^)]' + O(u-i6/0 UC 


as u —)■ oo. 


3. The asymptotic expansion of two integrals 


To evaluate the integrals appearing in (2.7) we employ the modification of Watson’s 
lemma given in Giver’s book [2, p. 337]. This is possible since Ai(x) possesses the 
asymptotic behaviour 

Ai(x) ~ exp[—|x^'^^] (x ^ +oo). 

Inversion of (2.3) with the help of Mathematica yields^ 


x = l + 2-i/=^(- 


-2/3^2 11(3 

To Too" 


823 • 2-^/3(/4 
12600 


+ 


(ICI<Co), 


^The function ^(x) is real for x > —1, being given by (—= |(arccos x — x+l — x^) in the interval 
— 1 < a; < 1. The finite radius of convergence results from the singularity in the inversion at x = —1, 
where C = 
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where Co = see [3, Eq. (12.10.41)]. Then we have 


<^(0 

OO 

= E 

m=0 

(ICI < Co), 

(3.1) 

where 

Oo — 2 ai — 

5 

25/3 

2i°/3 ^ 1548 


“ 35 ’ 

“ 225 ’ “ 67375 ’ " ' 



Now for a fixed positive integer M, and with t = 


oo 1 /*oo 

I </)(C) = XI 4 "Id / + i2M, 

0 ^0 Jo 


where 


M-l 


„C)0 

Rm = / </>m(C) Ai2(z.4/3^) ^ 

•^0 m=0 


It follows that 4>m{C) = O(C^) as C 0+ and it is known [1] that (p{() < 2“^/^ for 
C > 0. Hence we can find a constant C such that the inequality |C'm(C)I < holds 
when C ^ 0. Then 


\Rm\ <C r dC = 

Jo 

as z/ —)■ oo, upon use of the evaluation [3, Eq. (12.11.15)] 



t^Ai^{t)dt 


yjTiT{\m+l) 


{m = 0 , 1 , 2 ,...). 


Thus we obtain the asymptotic expansion 


(3.2) 


POO 

/ </>(C)Ai^(z^'‘/^C)'^C ~ 

Jo 


2 ^ a^ml I/-W3-4/3 


{v oo). 


(3.3) 


A similar reasoning can be applied to the second integral on the right-hand side in 
(2.7), when we observe that the function — 6o(C) ia (2.5) decreases monotonically from 
the value 9 • 2“^/^/140 to zero for C ^ 0 (we omit these details). With the help of 
Mathematica, the product (p{C)bo{C) has the series expansion 


mboic) = - E (ici < Co), 

m=0 


where the first few coefficients are 


/Jo 


,-1/3 


280 


/Ji 


iE2-2/3 g = 28687 
6300 ’ 2425500’ 


750979 

78828750 
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Then we find upon using an integration by parts, combined with (3.2) and the result 
Ai(0) = 3“^/^/r(|), that 


roo 

/ <('(C)feo(C) [Ai^ ~ 

Jo 


/?o(3u)-4/3 

r^(i) 




I3mm\ I/- W3-4/3 

12l”^+ir(im + |) 


as u —)■ oo. 


(3.4) 


4. The asymptotic calculation of Pn,tun for n —)• oo 


The function —ai(C) appearing in (2.4) is found to be monotonically decreasing for C ^ 0 
with the series expansion 


-ai(C) = 


249 


6849 I/O, 737 ^/oo 

2 "^^ C + 77^2-2/3(2-... (|C|<Co). 


28800 616000 65000 

Then the term multiplying Ai2(jy^/3() in (2.7) has the expansion 

-4<^(C)(1 + 1152/2(0) = x^2-2/3 - 25013 


576"^^^^" ^ 2400 1848000 

From (2.7), (3.3) and (3.4) we therefore find 


c + 


(ICI<Co). 


= ^ 1 + 


roo 

Jo 

29 \ 


2 aiu ^/3 2 q; 2 U ®/^ 603 ^ 

I ^ /ri / H~ //?-.—1/11 \ 


\/7rV 12zy2 2400zyWl2V6r(|) 12iV6r(^) 12W6r(^) 


+u-«/0l + 


1 \/3-4/3/lo ^ 2/3 iu-4/3 


+ 


121 / 2 ;vr 2 (|) ' ^/ri27/6r(|) 


+ O(z/-^®/0. 


Hence, from (1.2) and (2.6), we obtain 

2n+2^j^Qgn+2 


Pr 


n,tun — 


y/jT zy2n+5/3 


^ 1 29 \ f 6-2/3 jy-4/3 4 0 -l/ 3 j^- 8/3 

~ 24001/4 j I r 2 (l) ~ 3071^3 


525 r2(|) 525 r2(l) 

+ 0(U-1®/3)| 


+/,-8/3 1^1 , / A _ iA 

V 121/2^(280 r2(|) 6300 r2(l) 

upon use of the properties of the gamma function to simplify the coefficients. 
On noting that 


(4.1) 


2 ”+ 2 n!e ”+5 2^/3 5 


^^^2n+5/3 „l/3\^ 12^2 2881/^ 


23 


+ 0 {i 2 ®)^ (n —)■ oo). 
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we then finally obtain from (4.1) that 


Pr 


25/3 r 


n.tun — 1 /Q i 

n^/'^ I 


1 -- 


1 


g-2/3 j,-4/3 11 g_i/3^-8/3^ 


\ 167 6-2/3zy-4 _i6/3,l 


3i/Vvr2(i) 3071^3 600 r2(2) 


r^(3 


(4.2) 

as n —)• oo. Insertion of numerical values for the coefficients yields the alternative form 


P 


1 


n,tun - 1/3 


1 - 


1 


0.1339750 - 


0.0194484 

)r73 

0.0248598 


0.0174687 \ 

) 


+ 0(z/-i6/3)|^ i/ = \/2^^. 


The hrst two leading terms in this expansion can be seen to agree with the approximation 
stated in (1.3). 

We conclude by presenting the results of numerical calculations to illustrate the accu¬ 
racy of the asymptotic expression in (4.2). In Table 1 we show the values of the absolute 
relative error in the computation of Pn,tun using (4.2) for different n. 


Table 1: Values of the absolute relative error in the computation of Pn.tun from (4.2). 


n 

Pn^tun 

Error 

n 

Pfi^tun 

Error 

10 

0.0601438 

1.323 X 10"^ 

200 

0.0228302 

6.975 X 10"*^ 

20 

0.0483977 

2.528 X 10"® 

400 

0.0181454 

1.130 X 10"® 

50 

0.0360132 

2.534 X 10"'^ 

500 

0.0168499 

6.276 X 10"^® 

100 

0.0286973 

4.250 X 10"® 

800 

0.0144138 

1.815 X 10"^° 
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